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The effect of a strong magnetic field on the location of the critical end point (CEP) in the QCD 
phase diagram is discussed under different scenarios. In particular, we consider the contribution 
of the vector interaction and take into account the inverse magnetic catalysis obtained in lattice 
QCD calculations at zero chemical potential. The discussion is realized within the (2+1) Polyakov- 
Nambu-Jona-Lasinio model. It is shown that the vector interaction and the magnetic field have 
opposite competing effects, and that the winning effect depends strongly on the intensity of the 
magnetic field. The inverse magnetic catalysis at zero chemical potential has two distinct effects 
for magnetic fields above > 0.3 GeV^: it shifts the CEP to lower chemical potentials, hinders the 
increase of the CEP temperature and prevents a too large increase of the baryonic density at the 
CEP. Por fields eB <0.1 GeV^ the competing effects between the vector contribution and the 
magnetic field can move the CEP to regions of temperature and density in the phase diagram that 
could be more easily accessible to experiments. 


PACS numbers: 24.10.Jv, ll.lO.-z, 25.75.Nq / Keywords: PNJL, Polyakov loop, magnetic fields, critical 
end point 


I. INTRODUCTION 

The existence and the location of the critical end point 
(CEP) is a very timely topic for theoretical studies based 
on QCD with direct implications on the nature of the 
phase transition between the hadron gas and the quark 
gluon plasma. During the last decade, the investigation 
of the QCD phase diagram and the possible existence of 
the CEP has been attracting the attention of the physics 
community: by using lattice QCD (LQCD) simulations 
0,i, Dyson-Schwinger equations (sj and effective mod¬ 
els namely the Nambu-Jona-Lasinio (NJL) model Q, its 
recent extension up to eight quark terms (including ex¬ 
plicit chiral symmetry breaking interactions) [Sland the 
Polyakov-Nambu-Jona-Lasinio (PNJL) model ]^, there 
has been an effort to understand the nature of the phase 
transition and the existence of the CEP. 

From the experimental point of view the exis¬ 
tence/location of the CEP is the major goal of several 
programs namely the Beam Energy Scan (BES-I) pro¬ 
gram at RHIC which has been ongoing since 2010 look¬ 
ing for the experimental signatures of a first-order phase 
transition and the CEP by colliding Au ions at several 
energies Q . Recently, the results of the moments of net- 
charge multiplicity distributions were presented by STAR 
Collaboration Q. These measurements can provide rel¬ 
evant information on the freeze-out conditions and can 
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help to clarify the existence of the CEP. However, future 
measurements with high statistics data will be needed 
for a precise determination of the freeze-out conditions 
and to make definitive conclusions regarding the CEP . 
Also the dynamics associated with heavy-ion collisions, 
such as finite correlation length and freeze-out effects, 
should be considered in QCD calculations before dehni- 
tive conclusions about the CEP can be made Q . If there 
is a CEP with baryonic chemical potential, /i^, lower 
than 400 MeV, it is expected that the upcoming BES-II 
program can provide data on fluctuation and flow ob¬ 
servables which should yield quantitative evidence for its 
presence. Otherwise, late in the decade, the FAIR facility 
at GSI and NICA at JINR will extend the search of the 
CEP to even higher /i^ (for a review on the experimen¬ 
tal search of the CEP see [13 )■ Also the NA61/SHINE 
program at the CERN SPS aims the search of the CEP, 
and to investigate the properties of the onset of decon¬ 
finement through spectra, fluctuations and correlations 
analysis in light and heavy ion collisions [HI- 

There are several aspects that can influence the loca¬ 
tion of the CEP like the strangeness or isospin content of 
the in-medium or the presence of an external magnetic 
field [ 13 . Indeed, in Ref. [13, within the NJL model, it 
was verified that the size of the hrst order segment of the 
transition line expands with increasing B in such a way 
that the CEP becomes located at higher temperature and 
smaller chemical potential values. This was also verified 
by using the Ginzburg-Landau effective action formal¬ 
ism with the renormalized quark-meson model [l^ . The 
influence of strong external magnetic fields on the struc¬ 
ture of the QCD phase diagram is also very important 
because it has relevant consequences on measurements 
in heavy-ion collisions at very high energies [l^ . 
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At finite temperature, several LQCD studies have been 
performed to address the infiuence of the magnetic field 
over the deconfinement and chiral transition tempera¬ 
tures For a review in recent advances in the 

understanding of the phase diagram in the presence of 
strong magnetic fields at zero quark chemical potentials 
see 1^ . 

The inclusion of a magnetic field in the Lagrangian 
density of the NJL model and of the PNJL model allows 
describing the magnetic catalysis (MC) effect, i.e., the en¬ 
hancement of the quark condensate due to the magnetic 
field [23l - [^ , but does not describe the suppression of the 
quark condensate found in LQCD calculations at finite 
temperature and zero chemical potential which is due to 
the strong screening effect of the gluon interactions, the 
so-called inverse magnetic catalysis (IMC) (l6l - [l^ . In or¬ 
der to deal with this discrepancy, it was proposed that 
the model coupling, Gg, can be seen as proportional to 
the running coupling, a^, and consequently, a decreasing 
function of the magnetic field strength allowing one to 
include its effects {as{eB)). In fact, in the region of low 
momenta the strong screening effect of the gluon interac¬ 
tions weakens the interaction which leads to a decrease 
of the scalar coupling with the intensity of the magnetic 
field [131 ■ By using the SU(2) NJL model [1^ and the 
SU(3) NJL/PNJL models two ansatz were proposed 
that allow for the IMC. 

Other mechanisms that lead to the IMC can be found 
in the literature [l^, together with a model- 

independent physical explanation for the IMC while 
a review with analytical results for the NJL can be found 
in Ref. [s^. 

Another aspect that is relevant for the location of the 
CEP is the presence of the vector interaction which acts 
in the opposite way of the magnetic field [ 13 , In¬ 
deed, it is known that increasing the repulsive interac¬ 
tion strength in the quark matter phase diagram leads 
to a shrinking of the first-order transition region as the 
baryonic chemical potential increases (the CEP moves to 
larger and lower temperature T |39|). Furthermore, 
the increase of Gv can change the structure of the phase 
diagram by decreasing the possible quarkyonic phase M- 
It is also important to note that, as pointed out in 39|, 
there is no constraint for the choice of the coupling Gy 
at finite density; if we see Gy as induced in dense quark 
matter, the choice of Gy is not related with the vector 
meson properties in the vacuum but it can be related 
with in-medium modifications [43 |. 

The presence of a vector interaction also becomes rel¬ 
evant in reproducing some experimental results [43 | or 
compact star observations (see for instance d^), and so 
should also be taken into account in the computation of 
the equation of state (EOS) for magnetized quark mat¬ 
ter O, . A step toward this type of investigation has 
recently been taken in Ref. [33| where two flavor magne¬ 
tized quark matter in the presence of a repulsive vector 
coupling, described by the NJL model, has been consid¬ 
ered. 


In the present work we investigate the influence of the 
IMC and the vector interaction in the location of the 
CEP in magnetized matter using the (2-1-1) PNJL model. 
After the presentation of the model in Sec. [TTl several 
scenarios of interest will be explored: the influences of 
an external magnetic field and of the vector interaction 
in the location of the CEP when no IMC effects are taken 
into account (Sec. IIIII) : the influence of the IMC in the 
location of the CEP in the presence and in the absence 
of the vector interaction (Sec. IIVI) . Finally we draw our 
conclusions in Sec. 


II. MODEL AND FORMALISM 

The original PNJL Lagrangian is modified in order 
to take into account the presence of an external magnetic 
field and the vector interaction in ( 2 - 1 - 1 ) flavors: 

G = Q file] Q T Gsym T -^det T Gvec 

+ ( 1 ) 

where the quark sector is described by the SU(3) version 
of the NJL model which includes scalar-pseudoscalar and 
the’t Hooft six fermion interactions that models the axial 
1 /^( 1 ) symmetry breaking [13,1131, with £sym and £det 
given by 
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Aym = Gs'^ [{qKqf + i^lsKqf] , ( 2 ) 

a—0 

£det = -K {det [q{l + 75 )g] -k det [^(1 - 75 ) 9 ]} , (3) 
and a vector interaction given by [j^ ]. 
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£vec = -Gy [{qi^Xaq)'^ + (gT'^bSAa?)^] • (4) 

a—0 

Here, q = {u, d, s)'^ represents a quark field with three 
flavors, rhe = diagj(TO„, nid, to,,) is the corresponding 
(current) mass matrix, Aq = \/2j3I where I is the unit 
matrix in the three flavor space, and 0 < Aa < 8 de¬ 
note the Gell-Mann matrices. The coupling between the 
(electro)magnetic field B and quarks, and between the 
effective gluon field and quarks is implemented via the 
covariant derivative — iq/A^j^j^ — iA^ where 

qf represents the quark electric charge {qd = qs = 
-qu/2 = -e/3), and = d^A^^ - d^Af^ 

are used to account for the external magnetic field, 
and A^^{x) = gstrong A^{x)^ where is the SUc(3) 
gauge field. We consider a static and constant mag¬ 
netic field in the z direction, A^^ = 5^2XiB. In 
the Polyakov gauge and at finite temperature the spa¬ 
tial components of the gluon field are neglected: A^ = 
SqA^ = — 1^4 A^. The trace of the Polyakov line defined 
by $ = j^(('Pexpz dr Ai{x,T) ))^ is the Polyakov 
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loop which is the order parameter of the Z 3 symmet¬ 
ric/broken phase transition in pure gauge. 

The pure gauge sector is described by an effective po¬ 
tential U (<!),$; r) chosen in order to reproduce the re¬ 
sults obtained in lattice calculations 


a{T) 

T4 2 

-H 6 (r)ln [1 - -b 4(l>3 -b $3) - 3($$)2] , (5) 


where a (T) = oq -b Oi (^) -b 02 (^)^, b{T) = h ■ 
The standard choice of the parameters for the effective 
potential U is qq = 3.51, ai = —2.47, 02 = 15.2, and 63 = 
— 1.75. The parameter Tq is the critical temperature for 
the deconfinement phase transition within a pure gauge 
approach: it was fixed to a constant Tq = 270 MeV, 
according to lattice findings. 

As a regularization scheme, we use a sharp cutoff. A, in 
three-momentum space, only for the divergent ultraviolet 
sea quark integrals. For the parameters of the model 
we consider A = 602.3 MeV, niu = md = 5.5, MeV, 
i^ = 140.7 MeV, = 1.835 and KA^ = 12.36 as in 
. The thermodynamical potential for the three-flavor 
quark sector 12 is written as 


12(T, n) =Gs E {qfQff + 4:K (ququ) {qdqd) {qsqs) 

f—u,d,s 

-Gv Y. (9}9/)Ew(4>,$,r) 

f—u,d,s 

+ E + ^med + ^mag) (6) 

f—u,d,s 


where the vacuum 12™'^, the medium 12“®'^, and the mag¬ 
netic contributions 12™*^®, together with the quark con¬ 
densates {q/qf) have been evaluated with great detail in 
M- By minimizing the thermodynamical potential, Eq. 

with respect to the order parameters (q/q/), 4*, and 
<i>, we obtain the mean field equations. 

In the present study we consider the PNJL model with 
equal quark chemical potentials, /Tu = Mrf = Msj which 
corresponds to zero charge (or isospin) chemical potential 
and zero strangeness chemical potential (/rg = /is = 0). 

In [ 23 , it was shown that the running coupling de¬ 
creases with the magnetic field strength. Consequently, 
in the NJL model the coupling Gs, which can be seen 
as oc as, must decrease with an increasing magnetic 
field strength. Since there is no LQCD data available 
for as{eB), by using the NJL model it is possible to fit 
GsieB) in order to reproduce the pseudocritical chiral 
transition temperatures, T^{eB), obtained in LQCD cal¬ 
culations at /xs = 0 [I3- The resulting fit function that 
reproduces the T^{eB) is [ 2 ^ 


Gs(C) = G° 


/I + aC2 + K"\ 
U + cC^ + dCV ’ 


(7) 


with a = 0.0108805, b = -1.0133 x 10-^, c = 0.02228, 
and d = 1.84558 x 10“^, where ^ = bB/Aq^^y) ^-^d 
Aqcd = 300 MeV. 


By using GsieB) given by Eq. ([71) in the PNJL model 
both chiral and deconfinement transition temperatures 
decrease with increasing magnetic field strength due to 
the existing coupling between the Polyakov loop field and 
quarks within the PNJL model. Consequently, the cou¬ 
pling GsieB) affects not only the chiral transition but 
also the deconfinement transition 

Next, we will study the following scenarios for the effect 
of a static external magnetic field on the location of the 
CEP: 

1. Case I with no IMC effects and the usual Gg = G®: 

• Case lA, where we take Gy = 0 (Sec. IIII AI) : 

• Case IB, where we take Gy ^ 0, with Gy = 
aG° (Sec. Him . 

2. Case II with IMC effects, described by GsieB) 
given by Eq. ([7]): 

• Case IIA, where we take Gy = 0, (Sec. IIV AI) : 

• Case IIB, with Gy = aGsieB), meaning that 
the stronger the magnetic field the weaker the 
vector and scalar interactions, (Sec. IIVBI) : 

• Case lie, with a fixed Gy = aG° (Sec. IIVBI) . 


III. THE INFLUENCES OF AN EXTERNAL 
MAGNETIC FIELD AND OF THE VECTOR 
INTERACTION ON THE LOCATION OF THE 
CEP 


The influence of the repulsive vector coupling on mag¬ 
netized quark matter was investigated in by using 
the SU(2) version of the NJL model. Here, we will use 
the (2-bl) PNJL model. We start our study by setting 
Gy =0. In a second step we consider Gy ^ 0. 


A. Gv = 0 

The T — fjLB phase diagram obtained with Gy = 0 
(Case I A) is presented in the left panel of Fig. [T] and 
shows a trend very similar to that of the results previ¬ 
ously obtained for the NJL in [T^: as the intensity of the 
magnetic field increases, the temperature at which the 
CEP occurs increases monotonically (see Fig. [1] 

right panel) and the corresponding baryonic chemical po¬ 
tential (/r§^^) decreases until the critical value eB ~ 0.4 
GeV^ is reached; for stronger magnetic fields both 
and /r§^^ increase. In the middle panel of Fig. [1] the 
CEP is given in a T versus baryonic density, pb/P o, plot, 
and it can be seen that as the magnetic field increases 
from 0 to 1 GeV^, always increases. 

To understand these behaviors at finite density we 
start by considering the case at T = 0 where a first order 
phase transition takes place. In the left panel of Fig. [71 
we present the critical chemical potential, at which 
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FIG. 1. Location of the CEP on a diagram T vs the baryonic chemical potential (left), vs the baryonic density (middle), and 
vs magnetic field (right) for Case lA. 




eB(GeV") eB (GeV^) 

FIG. 2. The critical chemical potential at T = 0 MeV versus the magnetic field (left panel) and the u, d-quark constituent 
masses at the GEP as a function of the magnetic field intensity. 


the first order phase transition occurs. The pattern fol¬ 
lowed by at T = 0 in the PNJL model is similar, 
although for smaller values, to the one reported in [l^ 
at T = 1 MeV and also at higher temperatures: slow 
decrease for 0 < eB < 0.06 GeV^, faster decrease until 
0.12 — 0.18 GeV^, and monotonical increase afterwards. 
We verify a lowering of with B until eB = 0.25 
GeV^. The slow decrease in for increasing mag¬ 
netic field strength in the range 0 < eB < 0.08 GeV^ 
is followed by a faster decrease for 0.08 < eB < 0.25 
GeV^. Stronger field strengths result in a monotonically 
increasing This change in behavior corresponds to 

the point where just one Landau level (LL) is filled for 
each flavor in the partially chiral restored phase. Indeed, 
the stronger the magnetic field, the larger the spacing 
between the levels. 

At r = = 0 a stronger magnetic field results in an 


increase of the mass of the quarks (the increase is larger 
for Mu than Md due to the difference in electric charges). 
At finite density, however, starts to decrease with 
increasing magnetic fields, indicating an easier transition 
to the partially chiral restored phase [s^. This result 
was already seen in [I^. For eB above 0.25 GeV^, /r™* 
increases. 

Also noteworthy to point out is the existence of a range 
of magnetic fields, 0.083 < eB < 0.1 GeV^, where at least 
two first order phase transitions occur (see Fig. [21 left 
paneH), in accordance with what was found in the SU(2) 
[33, and SU(3) NJL models [s^. This cascade of 
transitions will result in the existence of multiple CEPs 


^ Around eB 0.085 GeV^ a small third phase transition (not 
visible on Fig. o can be found on a very small range. 
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FIG. 3. Effect of the strength of the vector interaction on the location of the CEP on a diagram T vs the baryonic chemical 
potential (left) and T vs the baryonic density (right). 


at finite temperature. The CEP on which we focus most 
of our attention in the present and next sections is the 
one that subsists to the highest temperature. 

As was discussed above, in the weak magnetic field 
regime an increasing magnetic field results in a smaller 
critical chemical potential for the first order transition at 
T = 0, even if the quarks’ masses have already started 
to increase. As this corresponds to a shift of the first 
order transition line toward a smaller chemical potential, 
the observed decrease in follows naturally. This 

effect is dominant over that of the increase of the quark 
masses at the CEP (both quark masses at the CEP in¬ 
crease with magnetic field strength for eB < 0.125 GeV^) 
which should hinder the first order partial chiral restora¬ 
tion (see Pig. [21 right panel). A similar behavior is also 
obtained within the NJL model used in 

Above a critical strength for the magnetic fields, eB > 
0.125 GeV^, there is a clear asymmetry in the CEP quark 
mass response to an increasing magnetic field strength: a 
strong decrease in Md as opposed to the smooth increase 
in Mu (due to the charge difference the d-quark coupling 
to the magnetic field is weaker). This behavior is accom¬ 
panied by an increase of the baryonic density at which 
the CEP occurs (Fig. [H right panel). 

For stronger magnetic fields {eB > 0.4 GeV^) both 
and increase. This can be understood as a 

result of a decreasing number of occupied LL due to the 
large intensity of the field and the greater difficulty in 
restoring chiral symmetry. 


B. Gv / 0 

The role of the vector interaction in the PNJL was 
studied in detail in [H, [s^. The main conclusion was 
that the CEP can be absent when the value of the cou¬ 
pling Gv is greater than a critical value G^'*. With 
the present parametrization this critical value is when 


a ~ 0.71, i.e, Gy'* w O.TIG® (see both panels of Fig. [3]). 
As the value of the coupling Gy is increased from 0 to 
Gy'* the first order phase transition is weakened and the 
CEP occurs at lower temperatures and larger chemical 
potentials but smaller densities. 

When the external magnetic field and the vector in¬ 
teraction are taken into account simultaneously, the sce¬ 
nario becomes more complex. In our discussion we will 
fix a = 0.25 (Gy = 0.25GO) and eB = 0.09 GeV^ to 
compare to cases with Gy = 0 and Gy = 0.25G° at zero 
magnetic field. The results are presented in Fig. 01 left 
panel. In the following we will show that for the light 
sector the (2-1-1) PNJL model is in agreement with the 
SU(2) NJL model (see [s^ for details). 

We start by exploring in detail the phase diagram in 
the presence of a magnetic field eB = 0.09 GeV^ at 
Gy = 0 where two first order phase transitions in cas¬ 
cade occur at T = 0. We identify a decrease of /ig'* for 
both transitions (see black and red full lines of Fig. 0] 
left panel) when compared with the transition at zero 
magnetic field (green full line in the same figure). 

At T = 0 these transitions occur for a given 
at which the gap equations have two stable solutions 
(M/ and M/^) leading to the coexistence of two dif¬ 
ferent quark masses at the same pressure and temper¬ 
ature. The first transition occurs at = 1061 MeV, 
being = 378.7(374.8) MeV and = 

166.0 (158.9) MeV, and the second transition occurs at 
HB = 1083 MeV, where = 149.8(141.4) MeV 

and = 49.4 (49.0) MeV. 

At this point some other relevant aspects of these re¬ 
sults are worth highlighting: 

i) At ^B = 0 and eB = 0.09 GeV^ we have the mag¬ 
netic catalysis effect as M„ and Md are higher than 
the respective vacuum values (M,™" = A/Y"'’ = 
367.7 MeV); 
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FIG. 4. QCD phase diagram in a T — /xb plot for Cases lA and IB, not including IMG effects, discussed in the text, left 
panel. The critical chemical potential as a function of the magnetic field intensity at T = 0, right panel. The green thick lines 
correspond to no magnetic field, and dashed (full) lines were obtained with (without) the vector interaction. 


ii) At the first phase transition the values of 
= 166.0 (158.9) MeV are still far from the 
respective current values (to„ = ttIii = 5.5 MeV), 
and a second transition is needed to bring 
to values closer to which is consistent with 

a region where the chiral symmetry is partially re¬ 
stored. 

hi) At the second phase transition = 

49.4 (49.0) MeV and these values are already 
smaller than = 52.2 MeV at eB = 0. 

iv) If we take the chiral limit for the light sector, 

= rrid = 0, the restoration of chiral symme¬ 
try, = 0, does not coincide with the first 

phase transition: at /i^ = 1012 MeV there is a jump 
in the quarks’ masses from = 366.5 (362.5) 

MeV to = 36.1 (33.5) MeV but only at 

fiB = 1017 MeV we have Af„ = Mii = 0. This is 
a direct manifestation of the condensate enhance¬ 
ment by magnetic catalysis. 

v) At finite temperature the second transition exists 
only at low temperatures and turns into a CEP at 
around T ~ 20 MeV (see red full line, left panel of 
Fig. 0]). Increasing the temperature, the remaining 
first order transition will subsist until 

MeV and = 856 MeV, closer to the CEP for 
eB = 0 (r™P = 156 MeV, = 873 MeV [H). 

When the vector interaction is also taken into account 
in the presence of a magnetic field, two competing effects 
come into play: on the one hand, an increase of Gy at 
eB = 0 weakens the first order phase transition which 
leads to the disappearance of the CEP in the fiB axis 
(this can be seen in the left panel of Fig. S] where the 
green dashed line shows how Gy affects the first order 
phase transition); on the other hand, an increase of the 


magnetic field at Gy = 0 has an opposite influence in 
the first order transition and CEP location, moving the 
transition line to smaller chemical potentials, at least un¬ 
til eB = 0.4 GeV^, as shown by the black line {eB = 0.09 
GeV^) in the left panel of Fig. [D Besides, as we already 
saw, the presence of a strong enough magnetic field can 
drive multiple CEPs due to the existence of several first 
order transitions at T = 0. 

Now, we discuss the case Gy = 0.25G° starting with 
T = 0 in the presence of an external magnetic field. As 
B increases, several first order transitions in cascade take 
place, similar to the Gy = 0 results. The critical chem¬ 
ical potential at which the transitions occur is shown in 
the right panel of Fig. S] (dashed lines): there is a range 
of magnetic fields, 0.03 ^ eB < 0.11 GeV^, where two 
transitions occur and a range, 0.07 < eS < 0.09 GeV^, 
where three transitions in cascade coexist. A first con¬ 
clusion about the combined effect of B and Gy is the 
appearance of intermediate transitions for a wider range 
of magnetic fields. 

Let us fix once again eB = 0.09 GeV^, a scenario 
where we have three phase transitions: compared with 
the eB = 0 result (green dashed line), two transitions oc¬ 
cur at lower hb (black and blue dashed lines in Fig. [U left 
panel) and a third one at higher ^b (red dashed line in 
the same panel). At finite temperatures the phase tran¬ 
sitions will give rise to three CEPs. The CEP with the 
larger temperature, T^ep _ ^qq ]y[gY = 1044 

MeV, (black dashed line in Fig. [U left panel), has a lower 
temperature and larger chemical potential than the CEP 
at eB = 0. 

In Fig. [5] left panel, the location of the CEP obtained 
with Gy = 0.25 Gs is plotted for different values of the 
magnetic field, the thin arrowed line shows the direc¬ 
tion of increasing fields. Going along the direction of 
increasing fields, there is a first region corresponding to 
the weaker magnetic fields B < 0.1 GeV^, where the 
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FIG. 5. The location of the CEP for Case IB, with the inclusion of the vector interaction and no IMC, the results for Case 
IB with a = 0.25 and 0.71 are compared with Case lA, with no vector interaction. The thin arrows indicate the direction of 
increasing magnetic field intensity. 


CEP temperature decreases and CEP chemical potential 
increases. Next, for 0.1 < eB < 0.4 GeV^ the CEP tem¬ 
perature increases and the CEP chemical potential de¬ 
creases, and finally for stronger fields the CEP chemical 
potential increases while the CEP temperature remains 
practically unchanged. 

For B < 0.1 GeV^, increases slightly and the 

respective baryonic density decreases, the magnetic field 
having an effect that adds to the one of Gy- In this 
region Gy is dominant and weakens the first order phase 
transition. This effect was clarified in [s^: although the 
effect of Gy goes always in the direction of reducing the 
density at the first order phase transition, the magnetic 
field does not present a monotonic behavior as seen for 
the scenario without the vector term. 

Above B = 0.1 GeV^, the effect of B on the GEP lo¬ 
cation is close to the one obtained for Gy = 0, except 
that for eB > 0.4 GeV^ the CEP temperature keeps in¬ 
creasing if Gy = 0, while for a finite Gy this trend is 
strongly reduced and the temperature does not change 
much; compare black and red curves of Fig. fright panel. 


C /~» /~»crit 

. Cry — Cry 

Finally, we investigate the case G™* w 0.71G°. The 
results are presented in the right panel of Fig. [5] (blue 
curves). The effect of the field is to restore a first order 
phase transition. 

For 0.01 < eB <0.1 GeV^ a very complex structure of 
multiple first order transitions appears at T = 0. These 
multiple transitions will be washed out with the increas¬ 
ing of the temperature until just one CEP remains. How¬ 
ever, more than one CEP can occur at very close tem¬ 
peratures like the scenario found in [l^. Above 0.1 GeV^ 
this complex structure disappears (in the right panel of 
Fig. Owe start to represent CEP for eB > 0.1 GeV^). 


For 0.1 < eB < 0.4 GeV^ we verify that the larger eB the 
larger until a maximum T ^ 145 MeV is reached. 

In the same range, the CEP chemical potential decreases 
slightly for eB < 0.3 GeV^, and increases for stronger 
fields. 

Increasing further the magnetic field, i.e. eB > 0.4 
GeV^, the increases while the temperature does 

not change much until eB = 1 GeV^, showing also sta¬ 
bilization of the GEP temperature even if weaker than 
the one obtained for Gy = 0.25 Gg. The number of oc¬ 
cupied LL becomes quite small. Indeed, taking eB = 0.4 
GeV^, the number of occupied LL at the GEP decreases 
with increasing Gy. The behavior obtained is the result 
of a clear competition between the magnetic field that 
disfavors chiral symmetry restoration and the vector in¬ 
teraction with an opposite effect. 


IV. THE INFLUENCE OF THE INVERSE 
MAGNETIC CATALYSIS IN THE LOCATION OF 
THE CRITICAL END POINT 

In this section we will investigate the influence, in the 
location of the CEP, of the IMC effect observed in LQCD 
calculations at zero chemical potential. First, the dis¬ 
cussion will not include the vector interaction, Case IIA 
(Sec. IIV A|) and next the IMC effects will be considered 
together with the vector interaction. Cases IIB and IIC 
(Sec. HVTO . 


A. Gv = Q 

The effect of the IMC on the CEP location excluding 
the vector interaction. Case IIA, is presented in Fig. [S] 
(red points) in the T — plane (left panel) and in the 
T — Pb/P o plane (middle panel), for different intensities 
















FIG. 6. Location of the CEP, versus (left panel) and versus (right panel), for different intensities 

of the magnetic field and excluding the vector interaction, without IMC effects Ga = (red curve) and with IMC effects 
Gs = Gs(eB) (black curve). 




(MeV) 


FIG. 7. Masses of the quarks as a function of ps for two intensities of the magnetic field: eB = 0.1 GeV^ (left) and : eB = 0.5 
GeV^ (right), at the respective 


of the magnetic field, and in the T — eB plane (right 
panel). For comparison we include in the same figure 
the CEP location without IMC effects, Case lA (black 
curve). We clearly observe a different behavior between 
these two scenarios: at i? = 0 both CEPs coincide but, 
already for small values of B, the IIA CEP is moved 
to lower temperatures and chemical potentials, keeping, 
however, a similar behavior to lA until eB ~ 0.3 GeV^. 
The large differences start for stronger magnetic fields: in 
Case IIA the position of the CEP oscillates between T ss 
169 and T « 177 MeV while the chemical potential takes 
increasingly smaller values; in Case lA both values of T 
and ps for the CEP increase (see black curve, left panel 
of Fig. [HI). In the middle panel of Fig. [Hlthe position of 
the CEP in the T — ps/po plane is presented. Comparing 
Cases lA and IIA, it is found that the IMC effect on 


the CEP results on its shift to smaller temperatures and 
densities especially for higher values of the magnetic field. 

The reason for these behaviors lies in the fact that the 
weakening of the coupling Ga{eB) will make the restora¬ 
tion of chiral symmetry easier. Increasing the magnetic 
field is not sufficient to counteract this effect, as can be 
seen in Eig. [3 where we plot the quark masses {Mu'- black 
line; M^: red line; Mg'- blue line) as a function of /is for 
the respective yC^p _ gq g^ _ q ^ QeV^. 

At eB = 0.1 GeV^, left panel, Gg is barely affected by 
the magnetic field when IMC effects are included and the 
values of the quark masses are very close to each other 
for both cases: in Case IIA the CEP occurs at smaller 
temperatures and at near, slightly higher, chemical po¬ 
tentials. When eB = 0.5 GeV^, right panel, the quark 
masses in Case lA have increased with respect to the 
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eB = 0 case (due to the MC effect), making the restora¬ 
tion of chiral symmetry more difficult to achieve. How¬ 
ever, when Gs = Gs{eB), Case IIA, the masses of the 
quarks are smaller than their eB = 0 value (due to IMG 
effect), leading to a faster restoration of chiral symmetry 
at small temperatures and chemical potentials. 

Eventually, with the increase of eB the CEP would 
move toward fiB = 0, and the deconfinement and chiral 
phase transitions would always be of first order. 


B. Gvt^O 

In this subsection we discuss the role of the vector 
interaction in the location of the CEP of magnetized 
quark matter taking into account explicitly the inverse 
magnetic catalysis at finite T through the renormaliza¬ 
tion of the coupling Gg due to the presence of the mag¬ 
netic field. We will consider two scenarios, with or with¬ 
out an explicit dependence of Gv on the magnetic field: 
Gv = aGsieB) (Case IIB), and Gv = ctG° (Case IIC). 
For the constant a we will take a general value a = 0.25 
(left panels of Fig. [5]), and the critical value discussed 
above, a = 0.71 (right panels of Fig. [5|) . 

We will first consider a = 0.25 [see panels (a), (b) and 
(c) of Fig. [S]. In these panels, the red curves corre¬ 
spond to Case IB with Gy = 0.25Gg, already presented 
in Fig. El and the black and the blue curves are for 
Gs = Gs{eB) and, respectively, for Gy = 0.25Gs{eB) 
(Case IIB) and Gy = 0.25Gg (Case IIC). When the IMC 
effects are included in the definition of the scalar cou¬ 
pling, and the vector term is taken into account, it is 
seen that for small values of B the CEP moves to lower 
values of T and slightly higher values oi (although a 
bit lower than in Case IIA), in both Cases IIB and IIC 
[see panels [5Kb) and [5Kc), respectively]. Although the 
coupling Gs{eB), and consequently also Gy in Case IIB, 
is slightly affected by the magnetic field, the overall bal¬ 
ance between the contributions of the attractive and the 
repulsive interactions is almost unchanged. 

For 0.09 < eB < O.SGeV^ the behavior for all three 
cases with Gy ^ 0, IB, IIB, and IIC, is very similar to 
the one found with Gy = 0 (Cases lA and IIA): the 
critical temperature increases [see panel [5Kb)] and the 
critical chemical potential decreases [see panel |5Kc)]. 

However, differences occur for stronger magnetic fields. 
When eB > 0.3GeV^, the CEP moves to smaller chemical 
potentials while the temperature does not change much in 
Case IIB. On the other hand, in Case IIC, with a coupling 
Gy that does not change with B, the vector contribution 
becomes more important than the catalysis effect of the 
magnetic field, since this last effect has been weakened by 
a weak coupling Gs{eB). Consequently, as soon as Gg is 
sufficiently weak, the effect of Gy is seen in the decrease 

of^CEP 

[see paneljSKb)], while practically does not 
change [see panel |5Kc)]. 

When a = 0.71 [see panels (d), (e) and (f) of Fig. [5] the 
CEP, for all cases with Gy ^ 0, shows a similar behavior 


as the held is increased from eB = 0, when pCEP _ q 
M eV and eB ~ 0.3 GeV^: the CEP moves to lower chemi¬ 
cal potentials [see paneljSKf)] and the critical temperature 
to larger values [see panel |5Ke)]. When the IMC effect is 
included. Cases IIB and IIC, the curves overlap due to 
a weakened coupling Gg{eB). Above eB ^ 0.3 GeV^, in 
Gase IIB the CEP occurs for smaller values of while 
2 ^cep (^Qgg change much, a behavior also occurring 
for a = 0.25. In Case IIC, with Gy = 0.71Gs, the trend 
is different from the corresponding one obtained with 
a = 0.25: the coupling Gg becomes sufficiently weak with 
respect to the magnitude of Gy and the effect of the vec¬ 
tor coupling is seen in the decrease of T^ep weak 

increase of For eB > 0.8 GeV^ the effect of the 

magnetic held clearly overlaps the effect of Gy, and the 
GEP goes to higher temperatures and chemical potentials 
as in Case IB without the IMC effect (Gy = 0.71G°). 

It is also interesting to study the effect of the mag¬ 
netic held on the CEP baryonic density. In Fig. |5[ the 
CEP location is plotted as a function of this quantity 
for different scenarios discussed in the present and previ¬ 
ous sections. It is seen that the CEP baryonic density is 
not much affected by the held, the only exception if the 
different behavior of Case IB for the stronger helds, the 
CEP temperature being the property that distinguishes 
the different scenarios. 


V. CONCLUSIONS 

In the present work we have discussed the possible con¬ 
sequences of the IMC effect on the location of the CEP. 
The discussion has been performed within the (2-1-1) 
PNJL model with the possible inclusion of the vector in¬ 
teraction. Within this model the IMC effect is described 
by considering a scalar coupling that weakens with the 
increase of the magnetic field as proposed in [2^ . The de¬ 
pendence of the coupling on the magnetic field has been 
fitted to the LQCD calculations of the transition temper¬ 
ature at zero chemical potential, which show a decreasing 
crossover temperature with an increase of the magnetic 
field for intensities below eB = 1 GeV^. 

The main conclusion of the present work is that the 
IMG effect will have noticeable effects on the location of 
the GEP. If the vector interaction is switched off, the CEP 
will occur at increasingly smaller chemical potentials and 
at a practically unchanged temperature if the field sat¬ 
isfies eB > 0.3 GeV^. This behavior is contrary to the 
findings of [l^ with SU(3) NJL with constant couplings, 
where it was shown that above 0.3 GeV^, both and 
increase. Also the baryonic density at the GEP is 
affected: including IMG effects it increases only 1/3 of 
the expected if IMG effects were not considered, making 
the GEP much more accessible in the laboratory. How¬ 
ever, for weaker fields, eB < 0.3 GeV^, both scenarios 
give similar results. 

If the vector interaction is included, we must consider 
two scenarios: a strong enough vector interaction will 
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FIG. 8. The location of the CEP: (a) and (d) T versus hb, for different intensities of the magnetic field; (b) and (e) T as a 
function of eS; (c) and (f) /is versus eB, for the vector interaction a = 0.25 [left, (a), (b) and (c) panels] and a = 0.71 [right, 
(d), (e) and (f) panels]. The red and bine curves were obtained with no IMG effects and correspond to Cases lA and IB. The 
IMG effects and vector interaction were included in the black curves (Gv is fixed) and magenta (Gv weakens with an increasing 
eB). For more details see the text at the end of Sec. IIVI 


turn the first order deconfinement/chiral transition into 
a crossover or, on the contrary, the vector interac¬ 
tion is not strong enough to wash out the CEP, but moves 
its location to smaller temperatures and baryonic densi¬ 
ties and to larger chemical potentials. For the first sce¬ 
nario we have confirmed within the (2-1-1) PNJL model 
the findings of [13 obtained with the SU(2) NJL model, 
showing that a strong magnetic field would transform 
the crossover into a first order phase transition. With re¬ 


spect to the second scenario, we have shown that for suffi¬ 
ciently small fields, the repulsive effect of the interaction 
is stronger than the MC effect originated by the mag¬ 
netic field, and the CEP location occurs at smaller tem¬ 
peratures and slightly larger chemical potentials. The 
decrease of the CEP temperature could be quite large, 
for a = 0.25: g^ffgj-g reduction above 60 MeV 

when eB goes from 0 to 0.09 GeV^. Both effects corre¬ 
sponding to the two scenarios move the CEP to regions 
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FIG. 9. The CEP location on the T-ps plane. Cases IIA, 
IIB and IIC with IMC effects are compared with Case IB 
without IMC effects. Case IIA is the only scenario excluding 
the vector interaction. 

of temperature and density in the phase diagram that 
could be more easily accessible to HIC. 

For larger helds, 0.09 GeV^ < < 0.3 GeV^ the mag¬ 

netic field wins and decreases while increases, 
showing a behavior similar to the corresponding one in 
the absence of a vector interaction. Above eB > 0.3 
GeV^ the CEP chemical potential increases but the CEP 
temperature keeps practically unchanged. 

The joint effect of the vector interaction and the IMC 
effect of the magnetic field depends strongly on the mag¬ 


nitude of the magnetic field and whether the vector cou¬ 
pling becomes weaker with the magnetic field. We have 
considered two scenarios: a constant Gy and a Gy that 
weakens when the magnetic field intensity increases. In 
the second case the location of the CEP for very strong 
magnetic fields is not much affected by the vector contri¬ 
bution, the magnetic field defining the structure of the 
phase transition. If, however, the vector coupling is not 
affected by the magnetic field, the weakening of the scalar 
coupling with the increase of the magnetic field inten¬ 
sity leads to the dominance of the vector contribution, 
which translates into a reduction of the CEP tempera¬ 
ture. An overall general conclusion is the reduction of 
the CEP temperature when the vector interaction is in¬ 
cluded. This effect will be quite strong if the vector cou¬ 
pling does not become weaker when the magnetic field 
intensity increases. 
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